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Unidimensional model of the ad-atom diffusion on a substrate submitted to a standing 
acoustic wave I. Derivation of the ad-atom motion equation 
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The effect of a standing acoustic wave on the diffusion of an ad-atom on a crystalline surface 
is theoretically studied. We used an unidimensional space model to study the ad-atom+substrate 
system. The dynamic equation of the ad-atom, a Generalized Langevin equation, is analytically 
derived from the full Hamiltonian of the ad-atom+substrate system submitted to the acoustic wave. 
A detailed analysis of each term of this equation, as well as of their properties, is presented. Special 
attention is devoted to the expression of the effective force induced by the wave on the ad-atom. It 
has essentially the same spatial and time dependences as its parent standing acoustic wave. 
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I. INTRODUCTION 

While the semi-conductors industry extensively uses 
the lithography process to stamp the micro-devices at 
the nanoscale, research centers and laboratories have in- 
vestigated the self-assembling properties of materials to 
avoid this expensive and time consuming process. Most 
strategies to self-assemble materials at the nanoscale, es- 
pecially during the atomic deposition process of semi- 
conductors benefit from the elastic properties or from 
the structure of the substrate: the Stranski-Krastanov 
growth mode relies on the competition between the sur- 
face and elastic energies to organize the 3D-growthji£ 
buried dislocations networks in the substrate induce a 
periodic strain field at the substrate surface that drives 
the diffusion of ad-atoms^ and finally the use of pat- 
terned substrates (vicinal surfaces, holes or mesas) can 
create some preferential nucleation sites £~— 

An alternative approach to self-assemble materials at 
the nano-scale, the dynamic substrate structuring ef- 
fect has been recently proposed^ At the macroscopic 
scale, a sand bunch on a drum membrane excited at 
one of its eigenfrequencies self-structures by accumu- 
lating around the nodes or anti-nodes displacements 
of the membrane^ Transposing this concept at the 
nanoscale, we investigate the diffusion of an ad-atom on 
a crystalline substrate submitted to a standing acous- 
tic wave (St AW) 2^ Molecular Dynamic simulations have 
evidenced that the StAW structures the diffusion of the 
ad-atom by encouraging its presence in the vicinity of 
the maximum displacements of the substrate^ The typ- 
ical and relevant StAW wavelengths vary from few to 
hundreds of nanometers. Experimentally, the produc- 
tion of standing surface acoustic waves of a few hun- 
dred nanometers to microns wave lengths are nowadays 
available through the use of interdigital transducer—^ 
or optically excited nanopatterned surfaces^ whereas 
one does not know yet how to efficiently generate smaller 
wavelengths (few to tens nanometers) phonons. 

In this study, we propose to analytically study the dif- 
fusion of a single ad-atom on a crystalline surface sub- 



FIG. 1. Schematic representation of the model under study. 
The ad-atom (red) and substrate atoms (blue) are charac- 
terized by their coordinates x and Xj (j G {— N..N}) in the 
reference frame of the center of mass of the substrate. Note 
that, for clarity reasons, the ad-atom is not reported on the 
same horizontal line as the substrate atoms, but the model is 
unidimensional. 



mitted to a StAW. The goal of this study is to estab- 
lish the formalism and the dynamic equation that de- 
scribes the diffusion of an ad-atom on a crystalline sub- 
strate submitted to a StAW. In section UIl a generalized 
Langevin equation governing the ad-atom diffusion on 
a one-dimensional substrate is analytically derived from 
the Hamiltonian of the system (ad-atom+substrate). 
Sects, [ml El El and EU detail the different terms in- 
volved in this generalized Langevin equation, as well as 
their properties. 



II. AD-ATOM MOTION EQUATION 

We consider the diffusion of an ad-atom on a crystalline 
substrate submitted to a StAW with a wave- vector in the 
x-direction. Since the ad-atom diffusion is expected to 
be mainly affected in the x direction, we specialize to a 
system with one degree of freedom. The extension to a 
2D system to model a more complex StAW system (for 
instance, two StAWs with wave vectors in the x and y di- 
rections form a square lattice of nodes and anti-nodes) is 
straightforward, though analytical calculations may be- 
come tedious. Fig. [1] reports a sketch of the model under 
study, x and x_ N , ..., x N respectively design the positions 
of the ad-atom and of the 2N + 1 substrate atoms in the 
reference frame of the center of mass of the substrate. 
Following the work of Zwanzigi^ and related works ^r— 



we start with the Hamiltonian of the isolated system (ad- 
atom+substrate) : 



H, 



2m 



-$(x,x_ 



N 2 



■> X n)+Y1 W J -+ V ^b(x_ N ,...,X N ) 



where m,p and rrij,Pj are respectively the masses and 
momenta of the ad- atom and of the substrate atoms, 
V au h(x_ N ,...,X N ) and $(x,X_ N , ~..,X N ) the potential 
enegies of the substrate-substrate and ad-atom-substrate 
interactions. At this point, the generation process of the 
StAW has not been yet introduced, this will be done later 
on. 

The motion of the substrate atoms will be described 
in the harmonic approximation^ with the associated 
phonons of eigenvibration frequencies w n , normal coordi- 
nates Q n and momenta II„. 



^ 2m 



■+V su b(x_ N ,....,X N ) - 2^2 [ n " n " +ulQ n Qn] , 

n 

(2) 
where over-bar quantities are complex conjugate quan- 
tities and where the potential origin has been fixed at 
the equilibrium positions, V su h(x°_ N , ...., x° N ) =0. In this 
equation and in all the following equations, unless other- 
wise stated, the summations over the substrate atoms j 
are from — N to N, and those over the normal modes n 
from — N to N excluding n — 0. Note that, within the 
harmonic approximation, for an isolated substrate there 
is no substrate dilation with temperature nor energy ex- 
changes between the phonons. 

The substrate atom displacements, Uj = Xj — x®, 
around their equilibrium positions, x^, are thus given by: 






k ^Q n , 



(3) 



where k n is the wave vector of the n normal mode. 
Eqs. © and © we have: 



In 



K—n Kn-i ^—n 



w n , Qn = Q~n, and II rl = II_ n . (4) 



*i?\X, X_ N , , 



From Eqs. @ and ([3]) and performing a developement 
of the potential <I> to first order in the Uj's: 

X N ) = ^(X'X^N 1 ■■■■T X N ) 

j 3 

= Mx) 

+ o H [Qn-$„(X) + Q n .y n {x)\5) 



where 



o 5$ 



(6) 



^(x) = J2^e^§-(x,x°_ N ,...., <). (7) 



The interaction of the ad-atom with the substrate has 
been separated in two contributions. $q(x), the first 
one, appears as an external static force field. It is due 
to the frozen equilibrated substrate inter-atomic periodic 
potential. The second one, represents the interaction of 
the ad-atom with the phonons Q n , i. e. with the moving 
substrate atoms around their equilibrium positions. 
Eq. (fTJ) hence writes 



ir 



H = -^- + <S>o{x) 
Zm 



- ^ [Qn-^n(x) + Q„.y n (x)] 

n 

- J2 [n„fi„ + UllQnQn] , 



(8) 



Note that the coupling between the substrate and the 
ad-atom is linear in the phonon variables and non-linear 
in the ad-atom variable, i.e. the reverse situation of the 
one studied by Cortes et al~Z. 

To model the presence of a StAW in Eq. ©, we add 
a forcing term with the same F amplitude on two spe- 
cific normal variables of opposite wave vectors Qn cx and 
Q„ CI (= Q-n cx ). However, since our model does not con- 
sider any dissipation of the substrate vibration modes, we 
slightly detune the forcing frequency f2 nox = <^ri cx + ^n BX 
from the eigenfrequency w„ ox to avoid any resonance and 
subsequent divergence of the amplitude of the mode Q Uox ■ 
These two modes will be equally excited and thus, from 
basic forced oscillation theory^ one expects a forced os- 
cillation substrate displacement field proportional to that 
of the parent standing wave: 

2F 
u(x, t) = - cos[^ Kcx £] cos(fc„ ox a; + 77), (9) 

where M is the mass of the oscillator, 77 a phase depend- 
ing on the initial conditions and with: 



A 2 = Ql cx - uj 2 ncx = K cx + Suj na J< 



(10) 



We thus consider the following Hamiltonian for the sys- 
tem (ad-atom+substrate submitted to a StAW): 



v 1 

H = 2m~ + $0 ^ + 2 51 [Qn-^n(x) + Q„.#„(a 

71 

+ ^^2[ttnttn+0jlQ n Q n ] 



-{Qn cx + Qn^FcOS^n^t] 



(11) 



Note that, in Eq. (jTlJ , the addition of the StAW term 
makes the Hamiltonian time-dependent. In addition, 
the work of the operator to induce the StAW (the last 
term of Eq. (ITT10 is not null on average and leads to a 
monotonous increase of the average energy of the system 
(ad-atom + substrate). This would be the case even 
taking into account all the nonlinear terms we have 
omitted in Eq. (jlip. We however assume that despite 



*;J ~T Oi,—j 



with Si j the Kronecker 



this monotonous increase of the energy, the temperature Where Aj j 

of the system remains constant, either by considering symbol 29 . In Eq. (|12b[) . — & n (x) is the force on the sub 

that the substrate is infinite and has the behavior of strate normal mode n, induced by the ad-atom at posi- 



a thermostat, or by considering that the system is not 
totally isolated and coupled to an external thermostat. 

The dynamic equations derived from Eq. (|lip readi^i 
dQn " (12a) 

„ - f n (x) + A„,„ ex Fcos[n„„i] (12b) 



dt 

dn n 


XI, 


dt 




dx 
~dl ~ 


P_ 

in 



(12c) 



dp _ _^£o/ s _ 1 \p 
dt ~ dx [X) 2 ^ 



,— — (x) +Qn—T-(x) 

dx dx 



(12d) 



tion x. Solving Eq. (|12a|) and Eq. (I12b[) between to and 
t, the normal substrate coordinates read: 



Qn(t) = Q„(io)cos[ W „(t-t )] + 5^sin[u; n (t-to)]- f ^n(x(t') / m[iUn{t t,)] dt'+ fA^Fcosin^t' f 11 ^ 1 *' )] dt\ 



U}„ 



LU,, 



to 



(13) 



I 

where Q n (to) and n n (i ) are fixed by the initial con- The left hand side term of Eq. (|17p is the usual inertial 
ditions. An integration of the second integral and an term. On the right hand side, we distinguish four terms, 



integration by parts of the first one gives: 
Q n {t) = C n (t ) cos[uj n (t - t )] + D n (t ) sin[w„(i - t )] 



Uos[u n {t-t>)]d X{tl) dJK {x{t , ))dtl 



to 



dt dx 



F 



A n : n cx -r2C0s[fi nox t], 



with: 

C„(t )= Qn(t Q )+L 



F 

A 2 



,,.,,, . -j-j cos[J7 n( , x to] 



A,(*o) 



Hn(*o) 



A r 



A 2 w 



From Eqs. ((4]) and (JTJ) we have: 

C„ = C_„ and I>„ = D_„ 



(14) 

*n(g(*o)) 

W 2 

(15a) 

■sin[f2„ ox i ]. 
(15b) 

(16) 



which are successively: 

• the force induced by the effective crystalline poten- 
tial $ c ff(x), defined by: 



Using Eq. (|12c[) . (|12d|) and (fT^jl . we derive the general- 
ized Langevin equation governing the ad-atom diffusion: 

m^J = -^|S(x) - J* j(x(t),x(t'), t ~ t')f t {t')dt' + m 



+FsAw(x,t). 



(17) 



$ cS {x) = ^o{x)--J2 — *n(a;)*n(a 



(18) 



The properties of this potential will be studied in 
Sect. IVH 



• The friction term - f* j(x(t), x(t'),t - t')^{t')dt' 
that depends on the ad-atom velocity and on the 
memory kernel j(x, x' , t — t'): 

1 {x^\t-t') = Y J C ° SM l' t ' )) d ^ix) d ^')- 

n " 

(19) 

The properties of 7(2;, x', t — t') will be studied in 

sect, inn 



• the stochastic forc o 17 i 18 £(£): 

n 

+D n (t )sm[u> n (t- t )} 



C n (to) COS[u n (t - to)] 
d^n 



dx 



(*(*))• (20) 



This term depends on the initial conditions and ad- 
atom position and is a quickly varying force gener- 



ated by the substrate. The properties of this force 
will be described in Sect. [V] 

• The last term FsAw(x,t) is the effective force due 
to the applied forcing term at fi„ ox , i.e. the force 
FsAw(x,t) induced by the St AW on the ad-atom 
through the substrate: 

F S Aw(x,t) = ^(^r t (x) + ^^(x)) cos[fi„.J]. 



dx 



dx 



This force will be detailed in Sect. Mil 



(21) 



The three first forces, crystalline, friction and stochas- 
tic, exist even in the absence of the StAW excitation. 
They are the usual forces describing the dynamics of the 
atoms in a crystalline material. 

We have chosen to keep in Fsaw(x, t) only the forced os- 
cillation term at the pulsation f2„ cx . All the other terms 
depending on F have been included in the stochastic force 
£(t). They correspond to the answers of the oscillators 
Q ncx and <3n ox to the initial conditions at i = to- Since 
the normal modes of the substrate are undamped, these 
last terms are periodic and do not cancel. For damped 
oscillators, the terms depending on F in the stochastic 
force would correspond to a transient regime and would 
thus cancel, contrary to the forced oscillation term at the 
pulsation f2 nex . 



III. THE STAW FORCE 

To derive the expression of the force FsAw(x,t) in- 
duced by the StAW, we need to explicit the expression of 

I 



$„ ex (a;) in Eq. (l2~Tj) . Since interaction potentials depend 
only on the relative position of the interacting particles, 
so do $o and ^> n . ^> n writes then: 



* n (x) = ^J=e^^(x-x°_ N ,....,x-xl) 

= a n (x)e ik «*, (22) 

with <x n (x) defined as: 

„0x <9$ 



^M = y J_ e -«n(x-»J) |i (a . _ x 

\/ m j V X 3 



(23) 



Note that for an infinite crystal the a n (x) functions have 
the lattice periodicity.^ . In addition, a n (x) — oc— n (x) so 
that introducing the real ot r n {x) = $t(a n (x)) and imagi- 
nary a % n {x) = ^s(a n (x)) parts of a n (x), we have: 

a r n (x) = a r _ n (x) and <(ar) = -at n (x), (24) 

which leads to: 



dx 

with 



+ ■ 



dx 



= 2[g n (x) cos(k n x) + h n (x) sin(fc n x)], (25) 



9-n 



da r n 
dx 



Kn&n 



-( fe «< + ^)- (26) 



where g n (x) and h n (x) have the lattice substrate period- 
icity. The FsAw{x,t) force writes then: 



FsAw(x,t) 



2F 

A2 



cos(fi„ ox £) [ffn cx (a:) cos(fc„ cx a;) + h n ^(x) sin(fc„ ox a;)] 



(27) 



The comparison of Eq. (|27|) to Eq. ^ shows that, as 
expected, the SAW force on the ad-atom, induced by the 
standing surface acoustic wave through the substrate, has 
the large scale spatial and time dependence of the cor- 
responding standing wave. This dependence at 27r/fc„ ox 
scale has been exhibited in Molecular Dynamic simula- 
tions 9 of ad-atom diffusing on a substrate submitted to 
a standing surface acoustic wave. However, at a finer 
scale, x smaller than the lattice parameter, this force ex- 
periences an amplitude and a phase modulation due to 
the presence of the crystalline potential through the func- 
tions g n (x) and h n (x). 

At this point, it is instructive to turn to a particular case 
by specifying the substrate and the interaction potential 
between the ad-atom and the substrate atoms, especially 
in order to get an explicit expression of the functions 
a n (x) and thus of g n {x) and h n (x). We assume that 



r 



the substrate atoms have the same mass M and that 
the ad-atom interacts with each substrate atom through 
an attracting pair potential V va i T {x — Xi) that cancels 
at infinity. We choose for Vp a i r an exponential curve of 
extension a (roughly the pair interaction range), i.e. a 
potential expression, that is physically meaningful and 
that allows the derivation of analytical calculations. 



<P(x,x°_ N ...x° N ) = y / V pail (x 



I> 



(28) 
Where Vo is the bonding energy. Note that minima of 4> 
correspond to atoms substrate positions. We have x^ = 
ja where a is the lattice spacing and j £ [-N, N]. a n (x) 



then writes: 



an (x) = ^=J2 e ' iknix ~ x 0j) jrM-> 



3 J 



o 

-JV : 



,*-<) 



to infinity (N — >• oo) in Eq. (|29|) . we obtain 

To 



a„(x) 



erx/M 



E 

j=m + l 



e -ik n (x-ja) e - 






0\ rfVnai 



pair 



f/.T 



(X - X°) 



E 



—ik n (x — ja) — — — ^- L 

e 'e " 



^LV^ e -ik n (x-ja) jL 



M 



dx 



(29) 



V, 



o 



uVM 



^ — ?fc Tl r+ — 



— ik n r— — 



-ik n a i i 



-ik n a- 



(30) 



(31) 



a n (a;) appears then as a function of r(x) only, which 
To take into account the cliscontinutics of the derivative writes: 



of Vpair at its minima, we define mo(x) and r{x) respec- 
tively, the quotient and the rest of the Euclidian divi- 
sion of x by a: x — m^a + r, with too G [-N, +N] and 
< r(x) < a. too(x) tells us in between which potential 
wells [toqcj, (mo + l)a] the ad-atom is and r(x) where it is 



= J^-o^ffl-c^)^ 



cr%/M cosh(f ) - cos(fc„a) 



(32) 



From this expression of a n we deduce the following 
exactly in between. Extending the size of the substrate expressions for \JJ n , g n and h n : 

I 



9n(r) 
h„(r) 



V e^cosh(£)-cosh(^) .. moa 



aVM cosh(f ) - cos(fc„a) 

V 



cr 2 VM[cosh(f ) - cos(fc„a)] 

Vo 

cr 2 VM[cosh(f ) - cos(fc n a)] 



T T — CI 

cos[(fc„(r — a)] sinh( — ) — cos(fc„r) sinh( N 

a a 



T T — d 

sin[(fc„(r — a)] sinh( — ) — sin(fc n r) sinh( ) 

a a 



(33) 
(34) 

(35) 



Note that, since g n (x) and h n (x) in Eq. (|27p have the of Vp a ir- A more symetric expression can be obtained 

lattice substrate periodicity, we have g n {x) = g n (moa + through the r = r' + a/2 translation, with now —a/2 < 

r ) = 9n{f) an d h n (x) = h n (moa + r) = h n (r). One r' < a/2 (r 1 = corresponds to the mid position between 

can easily verify that F$aw (Eq. (1271) ) is a continuous two successive potential wells, located at r' — ±o/2): 
function of x, despite the discontinuity of the derivative 



J 



Fsaw{x, r'(x),t) = F saw (r') cos(ft„ ox £) cos[fc„ ox (a; - r') + <po(r')] = F saw (r') cos(fi„ ox i) cos[fc ncx a: + <p(r')], (36) 



with 



F M 



2Fn 



fcn„^Q 



sinh — — cosh \- sin 



2o 



a 



k n „,Q- i 2 a • i 2 r 

cosn — sinn — 
2cr a 



2 «<n 



= 2F cos(^)sinh(|-) 

ta,n((po(r )) = tan — j*p— coth — - tanh — . 
2 2(7 a 



1 



1 + tan 



2 

2 ^cx^ 



1/2 



coth — sinh — 
2 2a J a 



1/2 



Fn 



2V F 



(37) 

(38) 
(39) 



A 2 cr VM[cosh(f ) - cos(fc„ cx a)] 

I 

where F saw (r') and <^(r') = fo(r J ) — k Uox r' are respec- spatial dependence of FsAw(x,r'(x),t). Eq. ([55]) again 
tively the amplitude and the phase of the large scale evidences the large scale spatial and time dependence of 



the SAW. This point is also evidenced by evaluating the at the midway position between two successive potential 
force at the substrate atoms positions, r' — ±a/2, and wells, r' = 0: 



Fsaw( x , r ' = ± a /2, t) = F sinh( — ) cos(fc„ ox :r) cos(f2„ ex £) 

a 

F S Aw(x,r' = 0,t) = 2F cos(^^)sinh( — )cos(fc„ ox o;)cos(n„ ex i) 

Z la 



(40) 
(41) 
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FIG. 2. Top and bottom: maximum force induced by the 
StAW (t = Of^ 2 -] in Eq. (O) and midle: interatomic po- 
tential (Eq. (|28|) ~). as a function of x/a, for fc„ cx a = 27r/15 
and two values of a/a: 0.5 (black) and 1.5 (red). Top and 
Bottom: envelop curve at the substrate atom positions (blue, 
Eq. (|4ip ) and midway in between ( magenta, Eq. (|40| l'). 



IV. THE MEMORY KERNEL 

Let's now study the memory kernel 7(3:, x' , t—t') of the 
friction force (Eq. (p~9|) ) that depends on the a n functions 
through ^ n (x) Eq. (|2"2|> : 



1 (x,x',t-t') = J2 



COS(u> n (t - t')) d* rl 

lo 2 „ dx 



W^(x') 



Note that this memory kernel depends on the ad-atom 
position so that the dissipation term in Eq. (I17|) is nonlin- 
ear in the ad-atom variables! 15 i An explicit expression 
of 7 is out of scope. However, since the a n functions are 
periodic functions of period the lattice parameter a, we 
can make an evaluation of the kernel without taking into 
account their spatial variations. They are then replaced 
in Eq. (|19p by their mean value over the period a. This is 
equivalent to take into account only the first term, ci„(0), 
of their Fourrier expansion: 

7 (x, A t-t') « £ ^(t-y-^ ^^ 



with 



(42) 



On(0) 



a n {x)dx. 



(43) 



Again using the particular inter-atomic potential 
(Eq. (|28l)) with Eq. ([31 ]) or (p | one gets: 



Fig. [2] reports both the maximum force (t = 0[ n ] in 

Eq. (f3"rj|0 induced by the StAW and the interatomic po- 
tential (Eq. (|28p) as a function of x/a for /c„ ox a = 27r/15 
and two values of a/a: 0.5 and 1.5. The large scale spatial 
dependence in cos(fc„ ox x) of the force Fsaw{x, r'(x),t) is 
clearly evidenced, whereas the finer scale, between two 
successive potential wells exhibits the sinus hyperbolic- 
based dependence of the force evidenced in Eq. (l37t . As 
a increases, the amplitude of the wells of the ad-atom- 
substrate potential (Eq. (|2"5]l) and the amplitude of the 
variations of the force at both the large scale 27r/fc„ ox and 
at the fine scale a decrease: indeed, if the interaction be- 
tween the ad-atom and the substrate is less pronounced, 
the force induced by the wave on the ad-atom well will 
also be reduced on both fine and large spatial scales. 



a„(0) 



2ik n V 



aaVM(kl + l/a 2 ) 
Within this approximation, the memory kernel reads 



(44) 



j(x-x',t-t') 



E 



a 2 M 



k^cos[uj n (t - t') 



l + kla 2 



ik ri (x — x') 



(45) 



In the same spirit, we will use the Debye model^ which 
is well adapted for simple monoatomic lattices at inter- 
mediate temperatures, to describe the phonon dispersion 
relation, u> n — c s k n , where c s is the speed of sound of 
the substrate; and change the discrete summation to an 



integral: 



W 2 a 2 
a 2 Mc 2 s 

kD k 2 cos[c s k(t - t')]e tk( - x - x "> 

■kn 



j(x-x',t-t') 



[1 + k 2 a 2 ] 2 



g(k)dk, (46) 



where kp — — is the Debye wave number and g(k) = 
L/(2ir) the density of states in the reciprocal space, with 
L = 2Na the size of the substrate. Moreover, since the 
function fc 2 /[l + k 2 a 2 ] 2 is a peaked function centered at 
k = of extension 1/cr, and considering that a is gener- 
ally larger than a, the limits of integration are extended 
to oo. An integration by parts leads to the calculation 
of Fourier transform of Lorentzians and to the following 
approximated 7 expression: 



j(x-x',t-t') 



LV 2 



2c 2 a 2 Ma 



H 



\x - x' + c s (t - 1') 



H 



\x-x' -c s (t-t')\ 



(47) 



with H(x) = (1 - x)e 



The expression of the memory Kernel in Eq. (|47|) is an 



even function of 



and t — t' '. The dependence on 



x — x 1 is a direct consequence of the elusion of the depen- 
dence of a n (x) on x (at the scale a) (see Sect. lVI[) . We do 
not find for j(x, x',t — t') a simple exponentially decreas- 
ing function of \t — t'\ as usually assumed in most text- 
books 26 . However, we emphasize that the 7 expression 
in Eq. (J47I) crucially depends on the interaction poten- 
tial chosen (Eq. ((28)) ) and that Eq. (|47|) provides a rather 
crude estimation of 7(0;, x', t — t'): we have ignored the 
dependence of *S> n on the length scale a and the extension 
of the integral Eq. (|4^|) to infinity is a rough assumption 
(a J a is not in general very large compared to 1). 
In addition, from Eq. (|4"7"1) the correlation time appears to 
be of the order of a/c s . Knowing that a is of the order of 
magnitude of the lattice paramater, this correlation time 
is of the order of the inverse of the Debye frequency. 



V. THE STOCHASTIC FORCE 

In this section, we describe the properties of £(£), 
the stochastic force (Eq. (f20|) ). Since this force de- 
pends on the ad-atom position through the coupling term 
^r-(cc(i)), it represents multiplicative fluctuations^ Us- 
ing Eqs. (I15a[) and (|15bp . it writes: 



m 



E 

n 

n«(«o) 



F *„(x(t )V 

Qn(t ) + A n> „ ox — COs[f2„ ox t ] H 5 ) COS[u n (t - to)] 



A, 



f a, 

A 2 " 



sm{n nc J Q ] sw.[w n (t - t )] 



d.r 



(x(t)). 



(48) 



This force partially results from the initial state of the 
substrate. In that sense, our system is completely de- 
terministic. However, we have considered a quadratic 
approximation in Eq. ([2]) and a linear development of 
$ in Eq. ([5])- In a real substrate, the non- linear terms 
can hold and/or exchange some energy with the normal 
substrate modes and in addition the substrate is never 
completely uncoupled to the experimental set-up. To 
take into account these exchanges of energy without ex- 
plicitly describing them, we characterize the state of the 
substrate {Q ,11) at to using a probability distribution 
p(Q(to),n(to)), where Q and II are vectors whose coor- 
dinates are the variables Q n and n„. We suppose that 
the StAW forcing terms in Eq. (|TT|) initially switched off 
are switched on at to: the Hamiltonian describing our 
system at t < to is thus given by Eq. © . 
Besides, if we want Eq. (JTTJ) to be regarded as a 
conventional Generalized Langevin equation, the quan- 
tity £(i) ought to have the properties that are ex- 
pected for Langevin noise. Especially, its average is 
expected to cancel with respect to the probability dis- 



I 

tribution p(Q(to), Il(to))21 In order to satisfy this last 
requirement, we choose the following expression for 
p(Q(t ),ri(to)): 



p(Q(t ),n(t )) = z- 



-i3H s 



(49) 



where f3 = l/(fcsT), ks the Boltzmann constant, T, the 
temperature of a surrounding thermostat that mimics the 
coupling of the system with the experimental set-up and 
H s given by: 



H S (Q,H) 



1 






[n„II„ + UjlQnQr, 



2^ 

n 

+ 2 E [Qn-*n(x{t )) + Qn-^n(x(t ))] 

n 

Flo 2 



[Qn«+Qn«] — rr-cospWo]- (50) 



A 2 



H s describes the coupling between the substrate and the 
ad-atom at position x(to) and contains a term derived 
from the StAW force to take into account the initial 
conditions imposed by the StAW on the Q n variables 
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at t = to. The Hamiltonian H s is hence different from 
the H one (Eq. flSJ)) of the system for t < t i.e. the 
probability distribution p(Q(t ),n(i )) corresponds to a 
non-equilibrium (macro-)state of the system described by 
H coupled to a thermostat at temperature T. We will 
now establish the properties of the fluctuating force £(£) 
for the probability distribution Eq. (|49|) . 

The examination of Eqs. (|4"51) and (|50l) reveals that the 
appropriate variables are: 

#n = Qn H — 5- + A„ jriGX —j cos[fi„ ox i ] ■ (51) 
<x>i A 2 

With these variables H s and £(£) write: 
1 ^ 



H S (Q,U) = 



2^ 



tu n R n R n 



w~ A^ 



c*o] 



(52) 



ew 



E>) 



fl n (to)cos[w„(t-t )] 

Ilnfa) _ . ^ ^ 



a.',, 



■sin[f2„ ox i " 



sm[w n (* - *o)] 



(53) 



From Eq. ([4T)[) and ([52")) , variables n„ and i?„ appear as 
complex variables with centered Gaussian distributions of 
variance /3 _1 . Note however, that since n„ = II_„ and 
R n = R- n , all these variables are not independent. One 
can easily re- write Eq. ([52")) using a set of 2N independent 
variables (R n ,H n ) with n > 0: 



H s (Q,n) = J2 



n>0 



uj n R n R n 



*, 



Ar. 



F 



cos[ft„ cx £ ] 



(54) 



So that, for any two variables X and Y G {uj n R n ,IL n } 
(n > 0) , their mean values (X) are and their covari- 
ances ([X - (X)][Y - (Y)]) are (2/f3)5 XY . 

From which we deduce the stochastic properties of £(£) 



<«*)) - |l 



* '™c 



(»(*)) + —J — (aw) 



da; dx 

sin[ri„ ox £ ] sin[w„ cx (t - t )] 



(55) 



C(t,t') 



(m- < m >} my < m >]> 

1 y^ cos[uj n (t — t')] 



P 



= - 7 (x(t),x(t'),t-t'). 



(x(t)) 



dx 



dx 



■(x(t'j) 



(56) 



We recover in this last equation the fluctuation- 
dissipation theorem: this result is especially independent 



of the precise expression of the potentials $ and V su h 
in Eq. ([T]) as soon as this later can be approximated by 
Eq. ©. The same result has been demonstrated in a 
general frame by Zwanzig^ The non-null value of (£(£)) 
is related to the time-depend Hamiltonian (Eq. (fTTj) ) and 
more precisely to the initial conditions that are imposed 
by abruptly switching on the St AW term at to. The 
Hamiltonian H s Eq. ([52")) actually takes into account the 
initial conditions imposed by the StAW on the Q n vari- 
ables but not on the n n variables. As a consequence, the 
non-null value of (£(£)) is directly correlated to the ini- 
tial conditions imposed on the n„ variables. To recover 
that the average value of the stochastic force cancels, we 
impose that r2„ cx i = 0[7r]: this corresponds to switching 
on the StAW force at an extremum of the force. 



VI. THE EFFECTIVE CRYSTALLINE 
POTENTIAL 

The effective crystalline potential $ c ff(x) reads: 

$ eS (x)=$o(x)+A$ (x), (57) 

with A$o(z) = -J2 Y^^n(x)^„(x). (58) 

n n 

Using Eq. ([22]). A$ (x) writes: 



A$o(a:) = - V] — ^a n (x)a„(x) 



— 1< 

n " 



(59) 



A$o(x) is then a periodic function of the lattice. A$ (x) 
physically corresponds to the modification of the poten- 
tial seen by the ad-atom induced by the auto-coherent 
interaction between the substrate atoms and the ad-atom 
at position x. Such interaction also appears in the mem- 
ory kernel. Actually, both terms A$o (x) and the memory 
kernel derive from the integration by parts of the third 
term of Eq. ([13) leading to Eq. ([T3|). The term A$ (x) 
derived from the third term of Eq. (|14[) . corresponds to 
the static and instantenous modification of the substrate 
variables due to the presence of the ad-atom at posi- 
tion x, while the memory kernel derived from the fourth 
term of Eq. (fl"4")) . corresponds to the retarded effects, i.e. 
how the past positions of the ad-atom influence the sub- 
strate positions at present. Both quantities A$o(a;) and 
j(x, x',t — t') can be related introducing the function 
@(x,x',t-t'): 

8(x, x\ t-t')=j: C0SK( * - 0) *n(aO ^ (*), 

n n 

dAa>0(x) =-i[e(a;,a;,0) + e(x,x,0)], 



d.x 
j(x,x',t — t') 



dx 



1 (x,x',t-t'). 



An explicit expression of the spatial dependence of 
A$ ( a; ) can be obtained using the particular inter-atomic 
potential (Eq. (|2"8)l). and the \l/„ expression of Eq. (|33l) . 
in Eq. (El 



A$ (z) = - 



l 2 



2cr 2 M 



cosh 2 (-) + cosh 2 ( ) 



a 
r , ,r — a, ■r—v 






-2cosh( — ) cosh( 

a a 



E 



[cosh(^) -cos(fe n a)]' 
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(60) 



cos(k n a) 



-^ u>% [cosh(f ) - cos(A;„a)]' 



where the two sums are only numerical factors indepen- 
dant of x. We recover in Eq. ([501) that A$o(a;) is a peri- 
odic function of the lattice. 



VII. CONCLUSION 

Wc have studied the diffusion of an ad-atom on a sub- 
strate submitted to a StAW. We found that the ad-atom 
motion is governed by a Generalized Langevin equation: 



d x (i$e: 

l ~dW = dx~ 

+F SAW {x,t) 



(*) + £(*)- / l(*,x',t 



t')f t {t')dt' 

(61) 



I 

We have characterized each of the terms involved in this 
equation and have given them their analytical expression 
and most of the time, an explicit expression. A key- 
result is the expression of the force Fsaw induced by 
the StAW as a function of x and t. Fsaw essentially 
varies as cos(fc„ cx a;) cos(f2 riox £) where fc„ ox and r2„ ox are 
the spatial and angular frequencies of the StAW. How- 
ever, a deeper analysis exhibits that this force also varies 
on the crystalline substrate lattice scale. The next paper 
of this series is devoted to the study of the solutions of 
the equation Eq. (|6T|) . 
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